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Abstract 

We consider the perturbation problem of a Mode III interfacial crack. The perturbation is of ge- 
ometrical type and can be both perturbation of the crack faces and perturbation of the interface, 
which can deviate from the initial straight line configuration. Asymptotic formulae are derived for 
the first-order perturbation of the stress intensity factor. It is shown that, due to the unsymmetrical 
nature of the problem, the Mode III skew-symmetric weight function derived in Piccolroaz et al. 
(2009) is essential for the derivation of the correct asymptotic formulae. 

To illustrate the method, we present the numerical results for different geometrical perturbations 
of a half-plane interfacial crack in an infinite bimaterial structure. Discussion on the extension of 
the method to finite bodies is also presented. 
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1 Introduction 



Modelling of interfacial cracks was addressed in the fundamental papers by Willis (1971) and Hutchinson 
et al. (1987). Oscillatory asymptotics of displacement and stress fields near the crack edges represent 
an important feature of interfacial cracks in vector problems of elasticity. Singular integral equations 
for the displacement jump across the crack take into account the full coupling between the normal and 
shear modes of loading. 

The notion of weight functions for cracks, as stress-intensity factors associated with a point force load, 
was introduced in Bueckner (1985, 1989). Another approach, viewing weight functions as special singular 
solutions of homogeneous problems for cracks, was used in Willis and Movchan (1995). Analysis of 
weight functions, represented by special singular solutions of the homogeneous problem for an interfacial 
crack was performed in Lazarus and Leblond (1998) and Antipov (1999). Further studies, involving 
both symmetric and skew-symmetric weight functions, are included in Bercial-Velez et al. (2005) and 
Piccolroaz et al. (2007). In particular, as the three-dimensional crack advances, the crack front deviates, 
and for the case of in-plane deviations the paper Piccolroaz et al. (2007) incorporates the first-order 
analysis of the weight functions based on the full solution of the corresponding matrix Wiener-Hopf 
equation. It has been shown by Piccolroaz et al. (2009) that the lack of symmetry inherent in the 
interfacial crack problem induces the appearance of a skew-symmetric component of the weight functions, 
so that, for instance, the skew-symmetric loads generate non-zero stress intensity factors even in the case 
of two dimensions (plane strain, plane stress or antiplane shear). 

We note that the problem of finding the variation of the stress intensity factors induced by a small 
geometrical perturbation of a plane crack placed at the interface between two dissimilar elastic materials 
requires the use of not only the symmetric weight functions but also of the skew-symmetric components. 
This statement is valid not only for a full vector problem, but it is also true for relatively simple anti- 
plane shear formulations. This is shown in the present paper with reference to the case of Mode III 
deformation. 

We consider both perturbation of the crack faces and of the interface, so that, in its perturbed state, 
the crack faces occupy the lines 

7± = {(xi,x 2 ) : £2 = ei/>±(xi), -oo < xi < 0}, 

whereas the interface occupies the line 

C e = {(xi,x 2 ) ■ x 2 = ecf)(xi),0 < xi < oo}. 

Here, we assume that the functions ip±, <fi and their derivatives vanish within a finite neighbourhood of 
the crack tip. Such perturbation is referred to as the regular perturbation of the boundary. 

Using the symmetric and skew-symmetric weight functions for the interfacial crack, we derive asymp- 
totic formulae for the perturbed stress intensity factor, which characterize the stress in the vicinity of 
the crack tip. 

The paper is organized as follows. The formulation of the problem and the governing equations 
are outlined in Section [2] Section [3] introduces the symmetric and skew-symmetric weight functions 
for the interfacial crack; the same section outlines the regular perturbation asymptotic procedure. An 
illustrative example and applications of the asymptotic formulae are included in Section 01 

2 Problem formulation 

We consider a twodimensional bimaterial structure made of two dissimilar materials, joined along a plane 
interface, see Fig.O The two materials are assumed to be linear elastic and isotropic, with elastic moduli 
denoted by /z± (shear modulus) and v± (Poisson's ratio). A plane crack is placed along the interface, 
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which is assumed to be perfect, so that the displacement and traction components are continuous across 
the interface £ e . 

The cracked body is loaded by out-of-plane tractions: P± = y u±i9u ± /9n ± along the outer boundary 
r±, and p± = [i±du ± / '8x2 along the crack faces 7^, where u denotes the displacement component along 
X3 . The loading is self-balanced in terms of both principal force and moment vectors 



/ P+dx + 


I P-dx — / p + dx -+■ 


I 




Jr. J~,' + 





p_dx = 0. 



(1) 



Suppose that the plane crack faces and the plane interface are slightly perturbed: the perturbation 
is denoted by X2 = eV>±(#i) with supports I± for the crack faces and X2 = e<p{xi) with support / for 
the interface, where e is a small parameter. We are interested in the calculation of the corresponding 
perturbation of the stress intensity factor Km. 




Figure 1 : Geometry of the problem 



The problem is then formulated in terms of the Laplace equation 

Au ± (x 1 ,x 2 ) = 0, (2) 



with boundary conditions: 



and interface conditions 



V±q^±=P± onr±, (3) 
0x2 

+ - du+ ' dvr ^ 
,J ' + ~d~ = tJ '~~d~ 01 



We assume that the loading p± along the crack faces 7j_ is separated from the supports I± . 
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3 Solution of the perturbation problem by means of symmetric 
and skew- symmetric weight functions 

The solution of the perturbed problem, u ± (xi, x 2l e), can be expanded for small e as follows 

u ± (x 1 ,x 2 ,e) = u k (x 1 ,x 2 ) + euf(xi 1 x 2 ) + 0(e 2 ), (6) 

where Uq (x\,x 2 ) is the solution of the unperturbed problem (e = 0) and uf(xi,x 2 ) is the first-order 
variation of the solution with respect to e. 

3.1 Model problem for Uq {x\,x<i) 

The model problem for u a {xi,x 2 ) is obtained by taking the leading order asymptotics of the perturbed 
problem as e = 0, so that u (xi,X2) must satisfy the Laplace equation ([2]) with the boundary and 
interfacial conditions ((3j) ((5^ where the boundaries 7j_ and £ e are replaced by 7° and £°, respectively, 
see Fig. 13.11 
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Figure 2: Model problem for Ug {xi,X2) 



The solution with locally bounded elastic energy and vanishing at infinity can be easily found by 
means of Mellin transform in polar coordinates. The Mellin transforms for the displacement vector and 
the stress tensor are defined as follows 



/ u(r,0)r s-1 dr, a=\ <r{r,9)r s dr, 
Jo Jo 



(7) 



and they are represented by analytic functions in the strips < Re(s) < 0.5 and —0.5 < Re(s) < 0.5, 
respectively. 

Correspondingly, the inverse transforms are 



-1 pLUl^-lOG -1 />UJ2^-iOO 

u{r,9) = — u( Sl 8)r- s ds, a(r,9) = — a(s,9)r-'- 1 d8, (8) 

where u>i,oj 2 lie in the respective intervals. 
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With taking into account all boundary conditions, one can obtain the solution in the form (see 
Appendix IA.ll for details) : 



u ± (s,( 



(p+ - p-)cos(s9) (p~P+ + H+P-) sin(s0) 



(/!+ + /i-)s sin(7rs) n±(/i + + /i_)scos(7rs) 



(9) 



Introducing the symmetric and skew-symmetric loading, [pj = p+ — p- and (p) = +p-)/2 
respectively, we obtain 



u ± fs, 



sin(s6 



-<P>(*) + 



p±s cos(ns) 
The stress components are given by 



cos(s9) 



— fi—) sin(s6 



cos(sO) 

— r^rwM 

COS(7TS) 

sin(s 



(/!+ + (-i-)s sin(7rs) 2/i±(/i + + (j,-)s cos(7rs) 
P±sm(s9) — cos(s( 



HM (io) 



C0S(7TS) 



(P>(* 



(/.*+ + /Lt_) sin(7rs) 2(/i+ + /x_) cos(7rs) 
H±cos(s6) (p + — p,_) sin(s#) 



+ //_) sin(7rs) 2(/i + + /i_) cos(7rs) 



IpIOO, 



(11) 

(12) 



This allows to obtain the two-terms asymptotics of tractions ahead of the crack tip and crack opening: 

(13) 



v27T V27T 



MM = 



(i+ + p- ( 2K 



m rl/2 _ 2Am r3/2 | +0(r5/2); 



/2tt 



3V2tt 



respectively, where 



(14) 

(15) 
(16) 



Am = -y- jf {<p)(r) + | WW} ^ 1/2 dr, 

= ^ ^ {<P>(r) + f MM} r- 3/2 dr, 
in which rj = (p,- — fx+)/(p+ + fi-) is the contrast parameter. 

3.2 Model problem for uf{xi,x 2 ) 

To obtain the model problem for uf{x\,X2) we will follow the procedure described in Movchan and 
Movchan (1995) and expand the term along 7± and ( € , and also along ( e . 
We obtain: 



9u± 




0"o 




9n ± 




9x 2 


■t'2 = 










dn 




9X2 


3J2 



9m± 



* 2 =o± I 9a;i 



9X2 



z 2 =0± 



9xi 



.''2=0 ± , 



i=0±, 



|+0(e 2 ), x!<0, (17) 
+ 0(e 2 ), xi>0, (18) 



and 



u± (xi , 0) + e (xi , 0) + </>( Xl ) 



9x9 



x 2 =0± 



+ 0(e 2 ), xi >0. 



(19) 
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Therefore the model problem for uf(xi,X2) (see Fig. 13. 2\ is formulated in terms of the Laplace 
equation (|2|) with boundary conditions 

duf 



along the boundary T±, 



M ± ^=0, (20) 



duf d , , . d ( , , , dufr\ 



along 7^., and the interfacial conditions 

»+ - „ r . 9(I1 ) := (g - g) - -to)g (l - Jt) . on C °, (22) 

Note that the right-hand sides in (|2ip. (j2"3")) can be interpreted as "effective" loading along the crack 
faces and along the interface, respectively, whereas (|22[) is a prescribed discontinuity of displacement 
along the interface. 

3.3 The reciprocity identity and evaluation of the stress intensity factor 

The integral representation of the Km stress intensity factor can be obtained by application of the Betti 
identity, in a way similar to Piccolroaz et al. (2007). The formula is extended here to include the case 
of imperfect interface with prescribed discontinuities of displacement u and traction a = fj.du/dx2- 
The reciprocity identity reads 

{lUl(x[ - sO^Xs!) + (U)(x[ - x 1 )laj(x 1 ) - (E)(xi - i,)H(n)}<iti = 0, (24) 
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where [[/J and (U) are the symmetric and skew-symmetric Mode III weight functions and (E) the 
corresponding traction (see Piccolroaz et al., 2009), 

— ^=x n 1 ^ 2 , xi > 0, 
lUj( Xl ) = { V2^ 1 (25) 

0, xi < 0, 

0, x% > 0, 

(E)(xx)= <J (l-i) M+M _ / ... (26) 

and (C/)(a:i) = 77/2[C/l(xi). 

Let us introduce the notations 



/(+)(x x ) = f{x 1 )H{x 1 ), f(-\ Xl ) = f( Xl )H(- Xl ), 
where H denotes the Heaviside function, so that 

/Ori) = /«(zi) + /«(zi). (27) 
The reciprocity identity (|2"3)) can be written as 

{ \U\ (x[ - xx ) (a) <+> (x x ) - <E) (xi - xx) H (") (xi) } da* = 

{[CKxi - xi)(a)(-'(^i) + (U)(x[ - a;i)H (->(*!)} <&i (28) 

) 

/>OG 

([/)(xi - x 1 )HW(x 1 )dx 1 + / (E)(xi - xx)H (+) (^i)^i- 



Note that (a)(~\ [[ct]|^ _ ^ are average and jump of the prescribed loading on the crack faces, whereas 

are the prescribed discontinuities of traction and displacement along the interface. 
The asymptotic procedure described in Piccolroaz et al. (2009) allows us to obtain the integral 
formula for the computatioix of Km as 

Km = -(l + i) Um / {lUj(x[ - xx)(a) ( - ) (x 1 ) + (U)(x[ - xx)H ( -)(xx)} dx x - 

poo poo 

(1 + i) Hm / (C/)(xi-xi)H (+) (^i)^i + (l + lim / (E)(xi - xi)[u] M (x 1 )dx 1 , (29) 



x^—tQjQ x'^ — >*0 _ 

which, upon substitution of ([25]) and ([^5]) , simplifies to 

Am = j=4 +tl - . / M (+) (^)^r 3/2 ^i- 

V27r(^ + + n_) Jo 

If {{a)^(x 1 )+ 1 llaf-\x 1 )}(-x 1 )- 1 / 2 dx 1 . (30) 

Note from (|30|) that the traction discontinuity prescribed along the interface, does not con- 

tribute to the stress intensity factor. This is because the weight function (U) is identically zero for 
xi < 0. 
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Now, we consider the asymptotic expansion of Km with respect to e as follows 

K m = K m + eK m + 0(e 2 ), 



(31) 



and make use of the integral formula (|30|) to compute both the leading term K® u and the first order 
variation K^ u . It is evident that the leading term K m corresponds to the unperturbed problem (e = 0) 
and was defined in (fT5")) . 

The second term K m corresponds to the model problem described in Sec. 13.21 and can be written as 



A III — TTT + TTT ; 



(32) 



where K m a ^ and K\ u ' denote the variations due to the perturbation of the crack faces and of the 
interface, respectively. 

Using (j3"D)) with a^' = /, and lu]^' = <?, we obtain 



K 



1(a) 
III 



^Xj-l)- 1 ^^ {</>(*!) + 

g(xi)x^~ 3/2 <ixi 



and 



A" 



,1(6) 



/27r(^ + + Jo 
where both jump and average notations are standard and 



(/)= *) 



9uo 
9xi 



1/1 = 



9uo 
<9xi 



, ^(xi) = -4>(xi) 



du 
dx 2 



in which the supports of the functions belong to the corresponding boundaries 7^. and £°. 
Integrating by parts the first integral (j33|) we get 



a iii — 



fiuoixx)— ((-x x ) 3/2 V>(xi)) 
in which, from the solution described in (j9l) . we can write for x\ < 0: 



dx\, 



1 / >io< 



f 2OO + 

where < <5 < 1/2, and for xi > 



1 =F r) cot its 77 tan 7rs 



2 s 



9x 2 



2TTlfl 



± J- 



COS 7TS 



(33) 
(34) 

(35) 



(36) 



(37) 



(38) 



4 An illustrative example 

To illustrate the perturbation method described above, we show in this section numerical results concern- 
ing the perturbation of both crack faces and interface geometry in the case of an interfacial half-plane 
crack in an infinite bimaterial structure, see Fig. |4] The crack faces are loaded by a "three-point" loading 
system consisting of: a point force F acting upon the upper crack face at a distance a behind the crack 
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tip and two point forces F/2 acting upon the lower crack face at a distance a — b and a + b behind the 
crack tip, see Fig. [4j In terms of the Dirac delta function #(■), the loading is then given by 



p+(xi) = F8(xi + a), 

F F 
P-{xi) = —S(xi + a + b) + —S(xi + a - b), 



or, in terms of symmetric and skew-symmetric parts, by 

(p)(*i) = ^S(xi +a) + -S(xi +a + b) + ^5( Xl + a - 6), 

IpKh) = <y(a?i + a) - -8{ Xl + a + b) - ^S( Xl +a-b). 
Here and in what follows we everywhere assume without lost of generality that F = 1. 

Xz" 



F/2 F/2 
"2b " 



2d +t 

A h 



^2d_ 



Material 1 



J\ \ A 
"2d " 



Material 2 



(39) 
(40) 



Figure 4: Illustrative example 



Applying the Mellin transform, we obtain 



{P)(s) = \a s + \{a + b) s + \{a-b) s , 
\p\(s)=a s -\{a+by-ha-by, 



(41) 



(42) 



so that, the zero-order stress intensity factor, corresponding to the unperturbed problem, is given by 

(43) 



^ III 



2 f'pa-^ + ^ia + by^ + ^ia-byVi 



4.1 Perturbation of crack faces 

Let us consider first the perturbation of the crack faces, defined by the two functions -0+ and -0-, which 
have the form 

i)±(xi) = ±-^r(^i + c± + d±f{xi + c± - d±) 2 , (44) 

and have supports in the intervals [— c± — d± , — c± + d±] . Note that the conditions c± +d± < a — b must 
be satisfied in order to have supports of the functions ip± separated from the loading. 



10 



From eq. (|3"6")l , we get 



K 



1(a) 
III 



1 

2^ 2 



c±-d± 



(^i)^-((-»i)- 3/ V±(^)) 



dxi, 



(45) 



where the term fJ.±u (xi) in the integrand can be computed using (|37|) . to obtain 



pi±wj(xi) 



IT?? 



, a 1 a + 6 I a-b. 

7 2 (log ) - -7 2 (log ) - -7 2 (log , 

— Xi 2 — Xi 2 — xi 



1+___ t n a 
— ^— -fi(log' 



-Xi' 



1 — n T a + 6 N l-?7., a — 6. . 



in which 



tan(-7rs) 



tanh(7r/j) 



1 

2?ri 



i r° 

h{0) = T -. / 

cotM A rtanhT^) A gi 



cos(/3t)dt, 



(46) 



(47) 



sin(j8i) ± tca&ifit) dt, (48) 



' —200+0 

where the sign depends on whether ±/3 > 0. 

The integrals 7i(/3) and 7 2 (/3) can be evaluated in closed form (see for instance Gradshteyn and 
Ryzhik, 2007): 

X - (49) 



h[fi) = -logcoth^ 

7T 4 



/a 08) 



/3 , 1 



2-7T 7T 



log 2 sinh ■ 



One can check that 



h((3) = — log^ + 0(/3 2 ), /3^0 
/_(,_) = l e -^/ 2 + 0(e-^), /3^±oo 
I 2 (/3)= .f + -log|/3|+0(/3 2 ), /3^0 

Z7T 7T 
ZTT 7T V / 



±00. 



(50) 

(51) 
(52) 
(53) 
(54) 



The estimations show that both functions 1\ and J 2 are singular for small argument values and also / 2 
tends to infinity for large positive values of its argument. 
The integral formula for K^f 1 is 



K 



1(a) 
III 



1±?7 



/2¥^ 2 



-c ± +a! ± 
-c±— d_ 



t n a \ 1 t n a + b \ 1 t n 

7 2 (log ) - -7 2 (log ) - -7 2 (log ) 

— Xi 2 — xi 2 — xi 



=F 



l + ?7,, a, 1 — 77 a + 6 , 1 — 77 a-b. 

—J-hOog ) + — -i/i(log ) + — -/i(log ) 

2 —xi 4 — Xi 4 — Xi 



^((-*i)~ 3/ V±(*i))<fci (55) 



Numerical results for different values of parameters are shown in Fig. __ We note that for the regular 

_____ l / 2 

perturbation geometry, the integrand in (|55[) may allow for a singularity not higher than x 1 ■ logxi 

for admissible values of the geometrical parameters a, b, c±,d±. 
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Figure 5: Numerical results for perturbation of the crack faces: (1 - green) rj = —0.99, (2 - orange) rj = —0.5, (3 - red) 
T) = 0, (4 - blue) rj = 0.5, (5 - black) rj = 0.99. 



Table 1: Parameters used for plots in Fig. \S\ 





Parameters defining the loading 


Parameters defining iji-y. 


Parameters defining t/>_ 


(a) 


a = l, < £>< 0.5 


A+ = 0.2, c+ = 0.4, d+ = 0.1 


A- -0.2, c_ = 0.4, d_ = 0.1 


(b) 


a = 1, < b < 0.5 


A+ = 0.2, c + = 0.4, (2+ =0.1 


A_ = 0.2, c_ = 0.4, d_ = 0.1 


(c) 


= 1, b = 0.2 


A + = 0.2, c + = 0.5, < d + < 0.3 


A_ = 0, c_ = 0.5, d_=0.1 


(d) 


a = l, 6 = 0.25 


< A + < 0.2, c+ = 0.4, d+ = 0.1 


A_ =0, c_ = 0.4, d_ = 0.1 


(e) 


= 1, b = 0.25 


< A + < 0.2, c + = 0.4, d + = 0.1 


A_ = 0.2, c_ = 0.4, d_ = 0.1 


(f) 


a = l, 6 = 0.25 


A + =0, c + = 0.4, cZ + = 0.1 


< A_ < 0.2, c_=0.4, d_=0.1 


(g) 


= 1, b = 0.2 


A+=0.2, 0.1<c+<0.7, d + =0.1 


A_ -0.2, c_=c+, d_=0.1 


GO 


= 1, b = 0.2 


A+=0.2, 0.1<c+<0.7, d + =0.1 


A_ = 0.2, c_=c+, d_=0.1 
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For parts (a) and (b) of this figure, we use skew-symmetric and symmetric perturbation of the crack 
faces, respectively. The geometrical parameters are shown in Table 1. The load is produced by point 
forces applied to the crack surfaces, as shown in Fig. @] A certain asymmetry is introduced in the 
applied load: one point force is applied on the upper crack face, whereas two smaller point forces at 
shifted positions are applied on the lower crack face; the total load is self-balanced, and the asymmetry 
is characterised by the ratio b/a representing the normalised distance between the two point forces on 
the lower crack face. Diagrams in Fig. [5^t,b show the first-order correction in the stress intensity factor 
for the cases of different asymmetry in the applied load and different values of the contrast parameter 
7/ = (//_ — + fi-). When the perturbation of the crack faces is skew-symmetric (see part (a)) 

the first-order variation in the stress intensity factor becomes zero for the case when the materials above 
and below the crack are equal and when b/a = 0, as expected. For the case of a bi-material interface, 
the sign of the perturbation depends on the material allocation and the stress intensity factor is shown 
to be sensitive to the asymmetry in the applied load. It is also visible on the diagrams of Figs. [5ji and 
[SJd that when the two forces are applied to the crack face on the stiffer half-plane (see the curve with 
the label 5) the sensitivity to the asymmetry in the applied load becomes insignificant. 

In Fig. [SJ;, we show the first-order correction for the stress intensity factor for the case when the 
lower face of the crack is not perturbed whereas the upper face has sustained a smooth perturbation 
with increasing support. It is shown that the increase in the size of the perturbation region leads to 
the increase in the stress intensity factor. It is also noted that the case of a soft upper half-plane, 
corresponding to the perturbed side of the crack, leads to the larger change in the stress intensity factor. 

In Figs. |3Ji and [S^, we show the cases when the lower face of the crack is not perturbed (as in 
Fig. [SJi) or perturbed with fixed amplitude (as in Fig. [5^), whereas the upper face is perturbed with 
the variable amplitude. It is shown that the correction for the stress intensity factor (initially zero for 
case[5ji, corresponding to the unperturbed geometry) grows linearly as the slope is different for different 
materials, following the same pattern as in the case [5h. Note that the values marked by the arrows in 
Fig. [5^, corresponding to A+/a = 0.2, are identical to the values marked by the arrows in Fig. 03, 
corresponding to b/a = 0.25, since the two cases correspond to the same geometrical configuration. 

Fig. [5f shows the case when the upper face of the crack is not perturbed whereas the lower face is 
perturbed with the variable amplitude, a situation opposite to the case[5ji. These computations confirm 
the conclusion that the variation of the stress intensity factor is stronger when the perturbation is applied 
to the softer material. 

In Fig. [Sg we show the results obtained for the skew-symmetric perturbation of the crack faces 
"moving" along the crack, so that the perturbed boundary becomes close to the region of applied tractions 
for c + /a 0.7~ and to the crack tip for c + /a —> 0.1 + . It is shown that the highest impact on the SIF 
appears when the perturbation is located near the crack tip, whereas a slight increase is observed when 
the perturbation appears near the acting load. It is also noted that for the homogeneous material this 
contribution is negligibly small. Fig. [5h corresponds to the case when the "moving" perturbation of the 
crack faces is symmetric. In this case, the influence of material parameters appears to be small (see also 
Fig. 03, where the influence of material parameters is enlarged by the refined scale) and the effect on the 
SIF increases considerably when the perturbation is approaching the crack tip, even for homogeneous 
material. 



4.2 Perturbation of the interface 

We consider now the perturbation of the interface, defined by the function 0, which has the form 

0(x 1 ) = ^ [ (x 1 +c + d) 2 {x 1 +c-d) 2 , (56) 
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and has support in the interval [c — d,c + d]. From eqs. and (pJTf . we obtain 



K 



i(b) _ FV 
in 



2V2tt 



-d 



(1 + ^(log A) + i-^J 3 (log ±±±) + i-^/ 3 (log — ) \ xf /2 <t>( Xl )d Xl , 



where 



w) = i stah| 



(58) 



7T sinh /3 

Numerical computations of the perturbation of the stress intensity factor for different values of 
parameters are shown in Fig. [6] The geometrical parameters characterising the perturbations are given 
in Table 2. 
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Fi gure 6; Numerical results for perturbation of the interface: (1 - green) r\ = — 0.99, (2 - orange) 77 = — 0.5, (3 - red) 
77 = 0, (4 - blue) 77 = 0.5, (5 - black) 77 = 0.99. 



Table 2: Parametersused for plots in Fig. [6] 





Parameters defining the loading 


Parameters defining <j> 


(a) 


o = l, 6 = 0.2 


A = 0.2, c = 0.5, 0<d<0.5 


(b) 


a=l, 6 = 0.2 


-0.2 < A < 0.2, c = 0.5, d = 0.25 


(c) 


o=l, 6 = 0.2 


A = 0.2, 0.25 <c< 0.75, d = 0.25 


(d) 


a=l, 0<6<1 


A = 0.2, c = 0.5, d = 0.25 



First, we note that the case of 77 = corresponds to zero perturbation of the stress intensity factor, 
since this situation corresponds to a straight crack in a homogeneous medium. In part (a) of the figure, 
we show the case when the amplitude of the perturbation is fixed, whereas the width of the perturbation 
region is subject to change. As expected, the increase in the size of the perturbation region leads to 
the increase of the magnitude of the stress intensity factor. If the tractions are applied to the crack 
faces at a finite distance from the crack tip then the asymmetry in the applied load does not give a 



14 



substantial influence to the magnitude of the perturbation of the stress intensity factor. On the part (b) 
of the figure, we look at the stress intensity factor as a function of the amplitude of the perturbation 
of the interface ahead of the crack. As expected, this dependence is linear, and the asymmetry of the 
applied load does not give much influence on the results of the computations. The diagram (c) of the 
figure corresponds to the case when the perturbation region of a fixed width and amplitude is displaced 
along the interface ahead of the crack. As expected, when the perturbation region moves towards the 
crack tip the magnitude of the perturbation of the stress intensity factor is increasing. Both cases, (a) 
and (c) show that the perturbation of the stress intensity factor is bounded, as predicted for the case of 
the regular perturbation of the interface. The last diagram (d) in Fig [5] illustrates the influence of the 
asymmetry of the applied load on the perturbation of the stress intensity factor. Namely, we consider the 
extreme situations when one of the forces applied to the lower face of the crack is allowed to approach the 
crack tip. In this case, Kjjj grows and correspondingly one can see the ratio of K)f) /K° IH decreasing. 
We also note that the curve with the label 5 corresponds to the stiff material in the lower half-plane, 
and hence the change in the stress intensity factor is visible only for small distance between the region 
of loading and the crack tip. 

5 Discussions and conclusions 

In this paper, we have illustrated the use of the skew-symmetric weight functions in evaluation of the 
stress intensity factors for regularly perturbed interfacial cracks in problems of anti-plane shear. The 
influence of asymmetry of the applied load as well as geometrical perturbations has been shown by 
the numerical simulations based on the explicit asymptotic formulae. Although the computations were 
presented for the case of a crack in an infinite two-phase plane, the symmetric and skew-symmetric 
weight functions for a finite domain can be constructed by a straightforward superposition of the infinite 
plane singular solution and an auxiliary solution with the finite energy constructed for the finite domain. 
The asymptotic algorithm for evaluation of the stress intensity factors, presented here, equally works for 
regularly perturbed domains of finite size. 
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A APPENDIX 



A.l Representation in terms of Mellin transform 

In this section we summarize the solution of the Mode III loading of the interfacial crack, by means of 
Mellin transform. 

The problem is formulated in terms of the Laplace equation 



d 2 ^ 9 2 u ± 



dx\ dx\ 



0. 



subject to the boundary conditions along the crack faces (x\ < 0) 



/'± 



dx 2 



x 2 =0 ± 



(A.l) 



(A.2) 



and the transmission conditions along the imperfect interface with presribed discontinuities of displace- 
ment and traction (x\ > 0) 



[u](a;i) = u + (xi,0 + ) - u (xi,0 ) = gi{x\) 

du 



ox 2 ox 2 



II- 



x 2 =0+ 



Ox? 



(A.3) 
(A.4) 



In polar coordinates the Mellin transforms for the displacement vector and the stress tensor with 
locally bounded elastic energy are defined as follows 



/•OO />OG 

u(s,6)= u(r,0)r s-1 dr, &(s,6)= a{r,6)r s dr, 
Jo Jo 



(A.5) 



and they are represented by analytic functions in the strips —$o < Re(s) < $oo and —70 < Re(s) < 700, 
respectively, where i?o 5 $oo > ($0 + 1?^ > 0), 7o,7oo > are constants related to the behaviour of the 
solution at the crack tip and at infinity, namely 



u{r,9) 



0{r°°) 
0(r-*« 



r -> 0, 



T — f 00, 

Correspondingly, the inverse transforms are 



tr(r,0) = 



O( r 7o-l) ) r ^Q ; 

O(r~ 7oo_1 ), t — y 00, 



u(r, 9) 



1 

2ni 



u(s,6)r s ds, cr(r, 



1 

2ttI 



0J2~\-ioc 



(A.6) 



(A.7) 



where — 1} < u>i < $oc an d —70 < w 2 < 700 • 

Applying the Mellin transform to the Laplace equation, 



Au± = - 



1 9u ± 9 2 u ± 1 9 2 u ± 



we obtain 



r dr dr 2 r 2 d9 2 
d 2 u ± (s, 



= 0, 



de 2 



+ s 2 u ± (s,6) = 0, 



The general solution of (|A.9[) is given by 

u ± (s,6») = A ± cos s0 + B ± sin s6 



(A.8) 
(A.9) 

(A.10) 
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In polar coordinates, the boundary and interface conditions take the form 

= p±{r), 



1 9m ± 



--±TT 



I(r) = u + (r,6 = 0+) - U -(r,9 = 0") = gi (r), 



and applying the Mellin transform 



1 du 



0=0+ 



00 



02 M, 



8=0- 



961 



P±(s), 



[u](a) = u + (s, = 0+)- fi-( S> = 0-)= 3i(s) 



m 



<9u 



0=0+ 



gi(r)r s 1 dr, 
g 2 (r)r s dr, 



so that 



The solution is given by 



fj,± (^sA^ 1 sin 7rs + si? 1 * 1 cos 7rs) = p± (s) 
A+-A- =g 1 (s), 
(ji + B+ - v-B-)s = g 2 (s). 

(P+ ~ P- i ^9l s sm + 02 COS 7Ts) COS S0 

+ /i-)s sin 7rs 

(— fA-P+ — /U+P- + H+H-9is sin7rs ± ^t±g2 cos7rs) sin s# 



(A.ll 

(A.12 
(A.13) 

(A.14) 

(A.15) 
(A.16) 

(A.17) 
(A.18) 
(A.19) 



(A.20) 

M±(/ i + + fJ--)s cosirs 

Introducing the symmetric and skew-symmetric loading, [pj = p+ — p— and (p) = (p+ +p-)/2 respec- 
tively, we obtain 



cos St 



— sin s6 



± 



(/t + + /Lt_)ssin?rs 2/i±([i + + /i_)scos7rs 
[i^ cos s9 [i + fi_ sin 7rs sin ,s0 



100 



sm sfc 



/i±S COS 7TS 



:<P>W 



(A.21 



/Lt+ + /i_ A ( ±(/ i + + /■*-) COS7TS 



5i 0) 



COS 7TS COS Sfc 



± 



+ /i-)s sin tts (/i++/i_)s 



In polar coordinates, the non-zero stress components are given by 



du ± 



9w ± 
<9r 



52(a) (A.22 



(A.23) 



and after Mellin transform, 



fi± sm sfc 



(/!+ — cos s6 



(/.t+ + /i_ ) Sin 7TS 2(/i + + /!_) COS7TS 



cos sfc 



(P>(«) 



(A.24) 



/i+/z_ssm; 



/Lt+ /i_ S Sm 7TS cos sfc 



+ /!-) COS 7TS 



0iO) 



/i-j- cos 7rs sin s0 fi± cos s9 



± 



+ (J,-) sin-KS + (j,-) 



92(a) (A.25) 
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-SjJL±U 



H± cos(s9) 



+ 



— sin(s6 



+ H-) sin(7rs) 2(/i + + cos(7rs) 



C0S(7TS) 



± 



fi+fj,-s cos s0 /i + /i_s sin 7rs sin s0 

/J+ + /i_ + //_ ) COS 7TS 



ffi(s) 



/x± cos 7rs cos s0 /i±sins0 



± 



+ |Lt_) sin7rs 



S 2 (s) (A.27) 



By means of the same procedure as in Piccolroaz et al. (2009), we obtain the asymptotics of stress 
and displacement near the crack tip as follows 



N 



<x(r, 0)=^2 Rcs[<t(.s, 6>), s = -kj2\r k l 2 - 1 + 0{r — ), 

k=l 
N 

u{r,6) = ^Rcs[m(s,6»),s = -k/2]r k/2 + O(r^). 



fc=0 



Since 



ct^(s,0) 



COS7TS 



{(P>W + fH( S )} 



1 1 
1 

/i_ J S COS 7TS I ' ' ' ' Z " * J COS 7T.S 

the two-terms asymptotics of tractions ahead of the crack tip and crack opening read 



^t + /!_S SUl 7TS _ /x± 

7 r <?l( s ) ± 7 \ 

[fl+ + /i_J COS7TS (^ + +/i_J 



52 (a) 



{<p>(«) + |ra(*)} 



H(r) = 



o"3e (r,0) 
2/^in ( \ 1 



2tt 



v 7 ^ \M+ ir- 
respectively, where 



■Kin -i/2 , ^ni 1/2 . -Bui 

—f=r 1 -\ -=r 1 + 

V2tt V2tt 

2 _2Am {_L + J_\ r 3/2 



.3/2 



0(r 5 / 2 ), 



(A.28) 
(A.29) 

(A.30) 
(A.31) 

(A.32) 



5\/27r \M+ M- 



in 



-/ {{P){r) + \ip\(r))r-^dr 



+ 



1 |U + /X_ 



.4 



in 
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in 



^ y o °°{(p)w+|ww}r- 3/2 *+ 

^°°{{p>(r) + |lp]W}r- B / 2 drH 



2S IXI f J_ + J_ ) r 5/ 2+0(r 7/ 2); (A 33) 



(A.34) 

(A.35) 
(A.36) 



2tt fJ,+ + fi- 



2tt /Lt + + /i- 
5 /i+M- 



2tt AH + M- Jo 



gx{r)r 3/2 dr, 
gi (r)r- 5/2 dr, 

o 

gi(r)r~ 7 / 2 dr. 
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